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■ Abstract. This paper studies a specific metric on plane curves that has the 

property of being isometric to classical manifold (sphere, complex projective, 
' Stiofel, Grassmann) modulo change of parametrization, each of these classical 

manifolds being associated to specific qualifications of the space of curves 
(closed-open, modulo rotation etc. . . ) Using these isometries, we are able 
' to explicitely describe the geodesies, first in the parametric case, then by 

modding out the paremetrization and considering horizontal vectors. We 
also compute the sectional curvature for these spaces, and show, in particular, 
1 that the space of closed curves modulo rotation and change of parameter has 

Q' positive curvature. Experimental results that explicitly compute minimizing 

' geodesies between two closed curves are finally provided 



(N 
> 



Introduction 



The definition and study of spaces of plane shapes has recently met a large 
amount of interest [2, 5, 7, 10, 18, 15], and has important applications, in object 
, recognition, for the analysis of shape databases or in medical imaging. The the- 

' oretical background involves the construction of infinite dimensional manifolds of 

, shapes [7, 15]. The Riemannian framework, in particular, is appealing, because it 

' provides shape spaces with a rich structure which is also useful for applications. 

, A general discussion of several classes of metrics that can be introduced for tliis 

' purpose can be found in [11]. 

. ^ , The present paper focuses on a particular Riemannian metric that has very 

' specific properties. This metric, which will be described in the next section, can 

^ I be seen as a limit case of one of the classes studied in [11], and would receive the 

label -ffi.oo in the nomenclature introduced therein. One of its surprising properties 
is that it can be characterized as a image of a Grassmann manifold by a suitably 
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chosen Riemannian submersion. A consequence of this is the possibility to derive 
exphcit geodesies in this shape space. 

A precursor of the -ffi^oo metric has been introduced in [18, 19] and studied in 

the context of open plane curves. It has also recently been used in [12]. Because 
the metric is placed on curves modulo changes of parametrization, the computation 
of geodesies naturally provides an elastic matching algorithm. 

The paper is organized as follows. We first provide the definitions and notation 

that wc will use for spaces of curves, the Hi ^c metric and the classical manifolds 
that will be shown to be isometric to it. We then study some local properties of the 
obtained manifold, discussing in particular its geodesies and sectional curvature. 
We finally provide experimental results for the numerical computation of geodesies 
and the solution of the related elastic matching problem. 

1. Spaces of Curves 

Throughout this paper, we will assume our plane curves are curves in the com- 
plex plane C. Then real inner products and 2x2 determinants of real 2-vectors 
are given by {x, y) = Bje{xy) and det(a;, y) = lm{xy). 

We first recall the notations for various spaces of plane curves which we will 
need, some of which were introduced in the previous paper [11]. For all questions 
about infinite dimensional analysis and differential geometry we refer to [9] . By 

Immop = lmm([0, 27r], C) 

we denote the space of C°° -immersions c : [0, 2tt] C. Here 'op' stands for open 
curve. Bi^op is the quotient of Immop by the group Difi^([0, 2tt]) of C°° increasing 
diffeomorphisms of [0,27r]. Next 

Immov('S'"'^, C), Immod('S'^, C) 

are the spaces of C°°-immersions c : S*^ C of even, respectively odd rotation 
degree. Here, is the unit circle in C, which will be identified in this paper to 
lR/(27rZ). Then Bi ^v and Si.od art; the quotients of Immov, respectively Immod 
by the group Difi^(5^) of C°° orientation preserving diffeomorphisms of S^. For 
example, Bi^od contains the simple closed plane curves, since they have index +1 or 
— 1 (depending on how they are oriented) . These arc the main focus of this study. 
To save us from enumerating special cases, we will often consider open curves as 
defined on but with a possible discontinuity at 0. We will also consider the 
quotients of these spaces by the group of translations, by the group of translations 
and rotations and the group of translations, rotations and scalings. 

Using the notation of [11], we can introduce the basic metric studied in this 
paper on these three spaces of immersions, but modulo translations, as follows. 
Identify Tc (Imm /transl) with the set of vector fields /i : 5^ — > C along c modulo 
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constant vector fields. Then we consider the limiting case of the scale invariant 
metric of Sobolev order 1 from [11], 4.8: 

(1) G,{h, h) = G™m,scal,l,oo(^^ h) = -^^ [ \D,h\^ .ds 

where, as in [11], ds = \co\d9 is arclength measure, Dg = Dg^c = \c0\~^de is the 
derivative with respect to arc length, and £{c) is the length of c. We also recall for 
later use the notation v = cg/\cg\ for the unit tangent vector, and, as multiplication 
by i is rotation by 90 degrees, n = i.v for the unit normal. Note that this metric 
is invariant with respect to reparametrizations of the curve c, hence it induces 
a metric which we also call G on the quotient spaces Bj^op,-Bi,ev and -Bj_od also 
modulo translations. 

The geodesic equation in all those metrics is a simple limiting case of those 
worked out in [11]. Suppose c{6,t) is a geodesic. Then: 

ctt = {{DsCt,v)DsCt - ^\DgCt\^v)) - {{DgCt,v)).Ct - ^{\DsCt\^) .D-^n.n) 

Here the bar indicates the average of the quantity over the curve c, i.e. {F) = 
J J Fds. Unfortunately, this case was not worked out in [11], hence we give the 
details of its derivation in Appendix I. The local existence and uniqueness of solu- 
tions to this equations can be proved easily, essentially because of the regularizing 
influence of the term D~^. This will also follow from the explicit representation 
of these geodesies to be given below, but because of its more general applicability, 
we give a direct proof in Appendix I. 

It is convenient to introduce the momentum u = — _Dg(ci) associated to a geo- 
desic. Using the momentum, the geodesic equation is easily rewritten in the more 
compact form: 

ut = -{u,DsCt)v- (^{DsCt,v) - {{Dsct,v))^u-^ ([D^ctp (|D,ct|2)) K(c).n. 

By tlK^ theory of Ricmannian submersions, geodesies on the quotient spaces Bi 
are nothing more than horizontal geodesies in Imm, that is geodesies which are 
perpendicular at one hence all points to the orbit of the group of reparametriza- 
tions. As is shown in [11], horizontality is equivalent to the condition u = a.n for 
some scalar function a{0,t). Substituting u = a.n and taking the n-component of 
the last equation, we find that horizontal geodesies are given by: 

at = -a ({Dgcu v) - {{DgCt, v))) + ^ (\DgCt\^ + {\DgCt\^)) . 

There are several conserved momenta along each geodesic t i— > c{6,t) of this 
metric (see [11], 4.8): The 'reparametrization' momentum is 

—^{cg,Dl^Ct)\cg\ 
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which vanishes along all horizontal geodesies. The translation momentum vanishes 
because the metric does not feel translation (constant vector fields along c). The 
angular momentum is 

Since the metric invariant under scalings, we also have the scaling momentum 

-i / {c,Dl^ct)ds = dtlog£{t). 

So we may equivalently consider either the quotient space Imm / translations or 
consider the section of the translation action {c e Imm : c(0) = 0}. In the 
same way, we may either pass modulo scalings or consider the section by fixing 
£{c) = 1, since the scaling momentum vanishes here. Finally, in some cases, we will 
pass modulo rotations. We could consider the section where angular momentum 
vanishes: but this latter is not especially simple 

2. The Basic Mapping for Parametrized Curves 

2.1. The bcisic mapping. We introduce the three function spaces: 
Vop = Vector space of all C°° mappings / : [0, 2'jt] — > M 

Vev = Vector space of all C°° mappings f -.S^ such that f{e + 27r) = /(6») 
Vod = Vector space of all C°° mappings f such that f{e + 27r) = -f{e) 

All three spaces have the weak inner product: 

/■27r 

Wff = / f{xfdx. 
Jo 

Given e, / from any of these spaces, the basic map is: 

$ : (e, /) ^ 0(9) = (1/2) / {e{x) + if{x)f dx. 

Jo 

The map c so defined carries [0, 27r] or to C. It need not be an immersion, 
however, because e and / might vanish simultaneously. Define 

Z{e,f) = {e:e{0)=f{0)=O}. 

Then we get three maps: 

{(e, f)€VxV: Z{e, /) = 0} ^ Imm,, for V = V^p, K^v, V^a. 

Looking separately at the three cases, define first the sphere S(yQp) to be the 
set of (e,/) e such that ||e||2 + y/f = 2. S°{V^p) is defined as the subset 
where Z{e, /) = 0. Then the magic of the map $ is shown by the following key 
fact [18]: 
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2.2. Theorem. $ defines a map: 

$ : S°(Vop) — > {c e Imiriop : i{c) = 1, c(0) = 0} = Imniop j {transl,scalings) 
which is an isometric 2-fold covering, using the natural metric on S and the metric 

(Jimm,scal,l,<x) Imniop. 



Proof. The mapping $ is surjective: Given c G Inrniop with c(0) = and £(c) ~ 1, 
we write c'{u) = r(u)e*"*^"-'. Then we may choose e{u) ~ ^2r{u) cos and 
/(w) = i/ 2r(u) sin Since 1 = 1(c) = f^^ |c'(w)| du — J^^ r{u) du we see that 
||e||2 + ||/||2 = /2^2r(w)(cos2(HM) + sin2(^iM)) du = 2. The only choice here is 
the sign of the square root, i.e. e, — /) = $(e, /), thus $ is 2:1. 

To see that $ is an isometry, let i>(e, f)=c = x + iy and (5c = (5a; + 1(5?/. Then 
the differential D^{e,f ) is given by 

(2) £>$(e, /) : (5e, Sf) ^ (5c(^) = y ' {Se + iSf) (e + if)d£. 

We have ds = {l/2)\e + if \^de. This implies first that £(c) = (||e||2 + ||/||2)/2 = 1 
as required. Then: 

_ 2(e + i/)(,5e + z(5/) 
Gc((5c, (5c) = i y |^.(<5c)|' ds = jj \6e + i5f\^d9 = \\{de, 5/)f . 



The dictionary between pairs (e, /) and immersions c connects many properties 
of each with those of the other. Curvature k works out especially nicely. We list 
here some of the connections: 



\ce\^U^' + f') 



ds 

de 

v = D(c)=^^±^ 
and if W${e, /) = efo — fee is the Wronskian, then: 

-e=[-^^)^ = ^^^^v, hence 
-2^^ for the curvature of c. 
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2.3. Geodesies leaving the space of immersions. Since geodesies on a sphere 
are always given by great circles, this theorem gives us the first case of explicit 
geodesies on spaces of curves in the metric of this paper. However, note that 
great circles in the open part S° are susceptible to crossing the 'bad' part S — S° 
somewhere. This occurs if and only if there exists 9 such that (e + if){9) and 
(^e + i5f{9)) have identical complex arguments modulo tt. So we find that our 
metric on Imm is incomplete. 

We can form a commutative diagram: 



where we have denoted the extended $ by <f>. For rather technical reasons $ is 
not surjective: there are pathological non-negative C°° functions which have no 
C°° square root, see [8], e.g. But what this diagram does do is give some space of 
maps to hold the extended geodesies. The example: 



- TT < X < TT, -1 < s < 1 and C = 2a^/3 + 2as^ 

and c(x, s) = {x^ — s^x + i3x^)/C — is/2 (suitably translated) 

is shown in figure 1. This is a geodesic in which all curves are immersions for 
s 7^ 0, but Cx{x, 0) has a double zero at x = tt. 

2.4. The basic mapping in the periodic case. Next, consider the periodic 
cases. Here we need the Stiefel manifolds: 

St(2, V) = {orthonormal pairs (e, f) £V xV},V = Vev or Kd 

and St" (2, V) the subset defined by the constraint Z{e,f) = 0. For later use, it 
is also convenient to note that St(2, V) = {A G L{C, V) : A^.A = Idc} where 
L(C, V) is the space of linear maps from C to V and {A'^v, w)c = {v, A'w)v- For 
(e, /) e S°(\4p), when is c = <I>(e, /) periodic? If and only if we have: 

(A) c' — i(e + i/)^ is periodic, so that e, / S Vcv or Vod, 



Condition (B) says that ||ejp = = 1 (since the sum is 2) and (e, /) = 

0, so that (e,/) e St°(2,Vev) or (e,/) € St"(2,Vod)- Recall that the index n 
of an immersed curve c is defined by considering log(c'). The log must satisfy 
log(c'(^ + 27r)) = log(c'(^)) + 27rn for some n and this is the index. So this index 




e{x, s) + if{x, s) = {x + is) / VC, where 
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Figure 1 . The generic way in which a family of open immersions crosses 
the hypersurfacc where Z 7^ 0. The parametrized straight hne in the middle 
of the family has velocity with a double zero at the bla<;k dot, hence is not 
an immersion. See text. 

is even or odd depending on whether the square root of c' is periodic or anti- 
periodic, that whether (e. /) G V^v or G Vod- So if $ is restricted to St°(2, T4v) or 
St°(2,V^(i) (and is still denoted $), it provides isometric 2-fold coverings 

$ : St°(2,Vov) — > {cG Imnicv(S'\C) : c(0) = 0,^(c) = 1} 

$ : St°(2,yod) {cG Immod(5'\C) : c(0) = (),£{c) = 1} 

All three of these maps $ can be modified so as to divide out by rotations. The 
mapping (e,/) 1-^ e*'^(e + if) produces a rotation of the immersed curve $(e,/) 
through an angle 2(p. The complex projective space CP(V^p) is S(V^p) divided by 
the action of rotations, and we denote by CP'^{V^p) the subset gotten by dividing 
S°(V^p) by rotations. The group generated by translations, rotations and scalings 
will be called the group of similitudes, abreviated as 'sim'. Then we get the variant: 

$ : CP°(V;2 ) Immop /(sim). 

Similarly, let Gr(2, V) be the Grassmannian of unoriented 2-dimensional sub- 
spaces of V and let Gr°(2, V) be the subset of those W with Z{W) = for F = V^y 
or = Vod- Then we have maps: 

$ : Gr°(2, Vov) — > Imm^^ /(sim) 

$ : Gr°(2, Kd) Immod /(sim) 
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For later use, we describe the tangent spaces of these spaces. The tangent space 
Tw Gr to Gr(2,V) at W € Gr(2,V') is naturally identified with L{W,W-^) and 
has the following norm, induced from that on V: 

\\vf = tv{v^ ov) = \\v{e)f + \\v{f)f 

for V £ Tw Gr and {e, /} an orthonormal basis of W. Similarly, T(^ej) St can 
be naturally identified with pairs {de,5f} in V such that {e,6e) = {f,5f) = 
(e, 5f) + (/, Se) = with norm: 

\\i6e,Sf)f = \\Sef + \\6ff. 

The same definition holds for T(gj)(S), this time with the constraint {e,Se) + 
{f,Sf)=0. 



3. The Basic Mapping for Shapes 



3.1. Dividing out by the group of reparametrizations. Let C°°'+([0, 27r]) 
be the group of increasing diffcomorphisms (p of [0, 2tt] (so (p{0) = 0, </?(27r) = 2tt) 
and let C°°'+(]R) be the group of increasing C°° diffeomorphisms : R ^ R such 
that ip{x + 27t) = (p{x) + 2-K for all x. Modulo the central subgroup of translations 
ip{x) — X + 2nn, the second group is just DifF"^ (5*-^ ) . For V = 14p,K;v, K)d let 
U{V) be the group of unitary maps on V given by: 



f ^ r = ■{! ° where 



(^G c°°^+([o,27r]) if y = yop 
(^e c~'+(R) if y = yev,Kd 



These are the reparametrization groups for our various spaces. The infinitesimal 
action of a vector field X on [0, 27r] or a periodic vector field X on is then 

(3) /^iXe./ + X/e. 



For all three sets of isometries we can now divide each side by the repara- 
metrization group U{y). For open curves, we get a diagram 

$ : S°(V;2p)/;7(v;p) Jul^ Bi^^^l ( transl,scalings ) 



* : CP°(V;2^)/f/(Kp) ^i?i,op/ (Sim) 
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and a similar one for closed curves of even and odd index where V = Vev, Vod and 
$ : St°(2, V)/UiV) B I (transl,scalings) 



$ : Gr°(2,y)/f/(y) ^Bj (sim) 

Here we have divided by isometries on both the left and right: by Uiy) or 
U (y) X on the left (where rotates the basis {e, /}) and by reparametrizations 
and rotations on the right. Thus $ is again an isoinctry if wc make both quotients 
into Riemannian submersions. This means we must identify the tangent spaces 
to the quotients with the horizontal subspaces of the tangent spaces in the larger 
space, i.e. those perpendicular to the orbits of the isometric group actions. For 
St, this means: 

3.2. Proposition. The tangent vector {6e, 6f} to St satisfies: 

{5e,e) = {Sf,f) = 0,{Se,f) + {Sf,e) = 0. 
It is horizontal for the rotation action if and only if: 

both Se, Sf are perpendicular to both e, f. 
It is horizontal for the reparametrization group if: 

We{e,Se) + We{f,Sf) = 
where W0{a,b) = a.be — b.ae is the Wronskian with respect to the parameter 6. 

Proof. Consider the action of rotations, which is one-dimensional, with orbits 
P e^^{e + if); the direction at (e,/) is chosen as (— /, e). So {Se,Sf) being 
horizontal at (e,/) for this action means that —{f,Se) + {e,Sf) = 0. This proves 
the first assertion. 

For the action of U{V), one has to note that horizontal vectors must satisfy 

{^Xe.e + X.ee, Se) + {\Xe.f + X.fe, Sf) = 

for any periodic vector field X on R, which yields the horizontality condition after 
integration by parts of the terms in X^i. □ 

Horizontality on the shape space side means (see [11]): 

3.3. Proposition, h S rcImm(S'^,C) is horizontal for the action of'DiS{S^) if 

and only if Dl{h) is normal to the curve, i.e. {v,Dl{h)) = 0. 

3.4. Proposition. For any smooth path c in Imm(S'^,M^) there exists a smooth 
path If in Diff(S'^) with ip{0, . ) = Id^i depending smoothly on c such that the path 
e given by e{t,9) = c{t,ip{t,6)) is horizontal: {Dl{et),ee) = 0. 
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This is a variant of [11, 4.6]. 

Proof. Writing instead ofD^ we note that D.oM^f) = \ cfo7 \ \Z \ = (-^«(/))° 
(p ior (p e Diff+(S'i). So we have L^^co^pif o ip) = {L^^cf) ° V- 

Let us write e = coip for e(i, 9) = c{t, ip{t, 9)), etc. We look for ip as the integral 
curve of a time dependent vector field ^(t, 6) on S^, given hj ipt = (p. We want 
the following expression to vanish: 

{Dlo^{dt{c o ^p)), de{c o ip)) = {Dl^^{ct o ip + {ce o ^p) ipt), {ce o cp) cpg) 

= {Dl{ct) o ip + DI(c0.£,) o p, C0 o p)(pg 

= {{{Dl {ct),ce) + {Dl {^.CB),cg))o^)^e. 

Using the time dependent vector field ^ = —j^D~'^{{Dl{ct),v)) and its flow ip 
achieves this. □ 



3.5. Bigger spaces. As we will sec below, we can describe geodesies in the 'clas- 
sical' spaces S, CP, St, Gr quite explicitly. By the above isometrics, this gives 
us the geodesies in the various spaces Imm, Bj. BUT, as we mentioned above for 
the space S, geodesies in the 'good' parts S°, CP°, St", Gr" do not stay there, 
but they cross the 'bad' part where Z{e,f) ^ 0. Now the basic mapping is still 
defined on the full sphere, projective space, Stiefel manifold or Grassmannian, giv- 
ing us some smooth mappings of [0, 27r] or to C, possibly modulo translations, 
rotations and/or scalings. 

But when we divide by U{Vop), a major problem arises. The orbits of U{Vop) 
acting on C°°{[0, 27r], C) are not closed, hence the topological quotient of the space 
C°°([0, 27r],C) by U{Vop) is not Hausdorff. This is shown by the following con- 
struction: 

(1) Start with a C°° non-decreasing map i/j from [0, 27r] to itself such that 
'tp{9) = TT for 9 in some interval /. 

(2) Let 'ilJn{9) = (1 — l/n).ip{9) + 9/n. The sequence {tpn} of diffeomorphisms 
of [0, 27r] converges to ip. 

(3) Then for any c G Immop, the maps c o are all in the orbit of c. But 
they converge to c o which is constant on the whole interval I, hence is 
not in the orbit. 

Thus, if we want sonic; Hausdorff space of curves which a) have singularities more 
complex than those of immersed curves and b) can hold the extensions of geodesies 
in some space Bi which come from the map $, we must divide C°°([0, 27t], C) by 
some equivalence relation larger than the group action by U{V). The simplest 
seems to be: first define a monotone relation R C [0, 2ti\ x [0, 27r] to be any closed 
subset such that pi{R) = P2{R) — [0, 2-k\ {pi and p2 being the projections on the 
axes) and for every pair of points (si,ti) € R and (52,^2) € R, either s\ < 82 and 
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ti < t2 or vice versa. Then /, g : [0, 27r] C are Frechet equivalent if there is a 
monotone relation R such that f{s) = g{t),y{s,t) G R. 

This is a good equivalence relation because if {/«}, {gn\ '■ [0, 2it] — > C are two 

sequences and lini„ /„ = /, lim„ g„ = g and fm9n are Frechet equivalent for all 
n, then f,g are Frechet equivalent. The essential point is that the set of non- 
empty closed subsets of a compact metric space X is compact in the Hausdorff 
topology (sec [1]). Thus if {i?„} are the monotone relations instantiating the 
equivalence of /„ and gn, a subsequence {Rnk} Hausdorff converges to some R C 
[0, 27r] X [0, 27r] and it is immediate that ii is a monotone relation making / and g 
Frechet equivalent. 

Define 

-Bbig.op = C°°([0, 27r], C) /Frechet equivalence, translations, scalings. 
Then we have a commutative diagram: 

$ : S°(V;2p) >Bi,op/ (transl,scalings ) 



$ : S(y„2p) ^ Bug,op 

Thus the whole of a geodesic which enters the 'bad' part of S(T4p) creates a path 
in Bug,o-p- Of course, the same construction works for closed curves also. We will 
see several examples in the next section. 



4. Construction of Geodesics 



4.1. Great circles in Spheres. The space S(y^) being the sphere of radius a/2 
on V^, its geodesics are the great circles. Thus, the geodesic distance between 
(e°, f) and (e\ f^) is given by s/2D with 

£> = arccos(((e0,ei) + (/°,r))/2) 

and the geodesic is given by 

^ sin((l-t)D) ^o sinfti^) ^! 
^ ' sinD sini? 

, S.n((l ,)D) s^^. 
sm U sm D 

The corresponding geodesic on Imm^p modulo translation and scaling is the time- 
indexed family of curves t ^ c{u, t) with 

dc/du = i(e(t) + imf = {e{tf - f{tf)/2 + ie{t)f{t) 
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The following notation will be used throughout this section; 

c°(u) = c(m,0), c^{u) = c{u,1) 

dc^/du = ro(u)e'"°("^ dc'^/du = ri(u)e*«'(") 

so that = y^2rj cos ^ and — y^2rj sin ^ for j — 0,1. Thus the distance D 
is 



-Dop(c'^, c^) = arccos / -^/roTicos du. 



i 



The metric on Imniop modulo rotations is 



-Don rot (c°,c^) =infarccos / y/ron cos — — —du 

"Jo 2 

/■27r 

Jo 



= arccos sup / \/ron" I cos cos — + sin sin — J du 

= arccos ^(^y V^r^cos ^ rfuj + y y^r^sin ^ dwj J 

The distance on op is the infimum of this expression over all changes of 
coordinate for c°. Assuming that c" and arc originally parametrized with l/27r 
times arc-length so that ro = l/27r,ri = l/27r, this is 

(4) Dop, difi(c°,c^) = arccos sup — / ^/^cos- 

and modulo rotations 

Dop, diff, rot(c°, c^) = arccossup (^(^^ y \/^cos ^ — —dd 



1 / r— . a — a ,„\-^ 



2n 



-de 



The supremum in both expressions is taken over all increasing bijections € 
C°°([0,27r],[0,27r]). 

To shorten these formulae, we will use the following notation. Define: 

1 r^'' r- a^o6-a° 



C_ (0) = — y^ y^^o cos ^ d0 

^-W = ^y["v/^sin^i^4^.. 

Then we have: 

-Dop, diff = inf arccos(C_ {(p)) and 
-Dop, diff, rot = inf arccos {^/{C^^yfT{S^^) . 
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4.2. Problems with the existence of geodesies. These expressions give very 
explicit descriptions of distance and geodesies. We have already noted, however, 
that, even if both (e°,/°) and (e^,/^) belong to S", the same property is not 
guaranteed at each point of the geodesic. e{a,t) = f{a,t) = happens for some t 
whenever (e°(a), /°(a)) and {e^ (a) , {a)) are collinear with opposite orientations. 
This is not likely to happen for geodesies joining nearby points. When it does 
happen, it is usually a stable phenomenon: for example, if the geodesic crosses 
S— S° transversally, as illustrated in figure 1, then this happens for all nearby 
geodesies too. Note that this means that the geodesic spray on Immop is not 
surjective. In fact, any geodesic on Immop comes from a great circle on S" and if 
it crosses S — S*^, it leaves Immop. 

When we pass to the quotient by reparametrizations, another question arises: 

does the inf over reparametrizations exist? or cquivalcntly is there is a horizontal 
geodesic joining any two open curves? In fact, there need not be any such geodesic 
even if you allow it to cross S — S*^. In general, to obtain a geodesic minimizing 
distance between 2 open curves, the curves themselves must be given paramctriza- 
tions with zero velocity somewhere, i.e. they may need to be lifted to points in 

s-so. 

This is best illustrated by the special case in which is the line segment from 

to 1, namely +i,f^ = l/v^i = 0- The curve c" can be arbitrary. Then the 
reparametrization </> which minimizes distance is the one which maximizes: 



i 



r- a° 
'(j)0 cos —d6. 
2 



This variational problem is easy to solve: the optimal (!> is given by: 

max I cos — , 



(u) = 27r y max ^cos — , 




Note that (p is not in general a diffeomorphism: it is constant on intervals where 
cos(a°/2) < 0. Its graph is a monotone relation in the sense of section 3.5. In fact, 
it's easy to see that monotone relations enjoy a certain compactness, so that the 
inf over reparametrizations is always achieved by a monotone relation. Assuming 
a° is represented by a continuous function for which — 27r < q^{u) < 27r, the result 
is that the places on the curve c° where |q!°(m)| > tt get squashed to points on the 
line segment. The result is that this limit geodesic is not actually a path in the 
space B of smooth curves. Figure 2 illustrates this effect. 

The general problem of maximizing the functional 

U{^) = ^1 V^cos " I d0 
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with respect to increasing functions <j> has been addressed in [17]. Existence of 
solutions can be shown in the class of monotone relations, or, equivalently func- 
tions (j) that take the form (/>(s) = /Lt([0, s)) for some positive measure /x on [0, 27r] 
with total mass less or equal to 2tt {(j)g being replaced by the Radon-Nykodim 
derivative of /i in the definition of U). The optimal (/> is a diffeomorphism as soon 
as cos((a^(w) — a^(v))/2) > whenever \u — v\ is smaller than a constant (which 
depends on a° and a^). More details can be found in [17]. 

It is easy to show that maximizing U is equivalent to maximizing 

[/+ W = ^ j^J v^max (^cos 0^ dO 

because one can always modify <f) on intervals on which cos((a^ o </> — a°)/2) < 
to ensure that (I)ed6 = 0. In [18], it is proposed to maximize 



o (h — a" 
cos 



d9. 



This corresponds to replacing the lift e{u) + ij{u) by a{u).{e{u) + ij{u)) where 
cr(u) = ± for all u, but this is beyond the scope of this paper. 




Figure 2. This is a geodesic of open curves running from the curve with 
tlic kink at the top left to the straight line on the bottom right. A blow up 
of the next to last curve is shown to reveal that the kink never goes away - it 
merely shrinks. Thus this geodesic is not continuous in the C-'^-topology on 
Bop ■ The straight line is parametrized so that it stops for a whole interval of 
time when it hits the middle point and thus it is C^-continuous in Immop. 
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4.3. Neretin geodesies on Gr{2,V). The integrated path-length distance and 
expHcit geodesies can be found in any Grassmannian using Jordan Angles [13] 
as follows: If Wo,M^i C V are two 2-diniensional subspaces the singular value 
decomposition of the orthogonal projection p of Wq to Wi gives ort ho normal bases 
{e°, f°} of Wo and {e\ of Wi such that p{e°) = Xe.e\ p{f) = A//i, e" ± 
and /° _L e^, where < A/, Ag < 1. Write Ag = cos('0e); A/ = cos{tpf ) then ^e, ^/ 
are the Jordan angles, < tpci'f < 7''/2. The global metric is given by: 



and the geodesic by: 



(5) 



W{t) 



sin((l - t).ipe)e° + sm{ttpe)e^ 



e{t) = 
f{t) = 



sin ' 
sin((l-t)%)/0 + sin(%)/i 
sin tjif 



Wc apply this now in order to compute the distance between the curves in 
the two spaces Immcv /(sim) and Imniod /(sim), as well as in the unparametrized 
quotients i3i,ev/(sim) and Bi^o^/ We write as above doc^ = ro(6')e*" and 
agc^ =ri(6')e'«'W. We put 



/2rosin ^ 

1 

'2ri sin ■ 



= \/2ri cos ^ ^ — v — i 2 

thus lifting these curves to 2-planes in the Grassmannian. The 2x2 matrix of the 
orthogonal projection from the space {e°, /°} to {e^, /^} in these bases is: 

7gi 2Vr°.rK cos ^ cos ^ '^6' /^i 2Vr0.ri. cos sin ^ 

Jgi 2\/^^. sin ^ cos 4 d6l 2\/^^. sin ^ sin ^ dO/ 
It will be convenient to use the notations: 

cos d0 



S± := 



■-L 

:= / VrO.ri 



sm 



i(M(cO,ci)n)TM(cO,ci)22) 
i(M(c°,ci)2i±M(c°,ci)i2) 



We have to diagonalizc this matrix by rotating the curve c" by a constant angle 
/S", i.e., the basis {e°,/"} by the angle and similarly the curve by a 

constant angle /3^. So we have to replace a" by a° — /3° and by — in such 
a way that 



(6) 



5±.cos^-C±.sin^ 



dO (for both signs) 
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Thus 

/3o ± /3i = 2 arctan (5±/C±) . 

In the newly aUgned bases, the diagonal elements of M(c°, c^) will be the cosines 
of the Jordan angles. But even without preliminary diagonalization, the following 
lemma gives you a formula for them; 

Lemma, /f M = ^ ^ and C± = ^{a^ d), S± = |(c± 6), then the singular 

values of M are: 

The proof is straightforward. This gives the formula 



(7) i?od,rot(c°, c^f = arccos^ J Si ^-C\^- J Si + Cl 



+ arccos^ [^^ Sl + Ci - ^ Si + ^ 
This is the distance in the space Immod(<S'^, C)/(transl, rot., scalings). 

4.4. Horizontal Neretin distances. If we want the distance in the quotient 

space i3i^od/(transl, rot., scalings) by the group Diff(5^) wc have to take the infi- 
mum of (7) over all reparametrizations. To simplify the formulas that follow, we 
can assume that the initial curves (P,c^ are parametrized by arc length so that 
r" = = 1/2-K. Then consider a reparametrization <j) € Diff(S'^) of one of the two 

curves, say cP o 0: 

(8) Aim,diff(c°, c^f = inf (arccos2(Ae(c° o </., c^)) + arccos^ (A/ (c° o </), c^))) 
where now 

Ae(cO o 0, ci) = ^Sl{^) + Cl{4>) + ^Sl{(t>) + C2 (</,) 
A/(c° o cj>, c') = ^Sl{cl^)+Cl{<l>) - ^J^W) + Cl{^) 

5±(</.) := ^ j^^ v^sin i^!£^ de, 

C±{ct>) ■.= ^ I V^ecos (2!^f±£l de. 

To describe the inf, we can use the fact that geodesies on the space of curves 
are the horizontal geodesies in the space of immersions. Consider the geodesic 
t ^ {e{t), f{t)} in Gr(2, V) described in (5), for 

e^ = ^cosi^^ e^ = ^cosi^, 

f° = ^f^sm ("°°t^-^° /I = ^ si" 



Sin ■ 
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where the rotations f3^ and (3^ must be computed from o (j) and . Note that 



2 



2^(/)/ .(aeO</,).cos- 

2 

If the Jordan angles are ■(/'e and ■0/, then the tangent vector to the geodesic t ^ 
W{t) at t = is described by 

e,(0) = atloe = (e^ - cos^e-e") , /t(0) = 5*|o/ = (/^ - cos^/./°) 

By 3.2 the geodesic is perpendicular to all Diff (S'^)-orbits if and only if the sum 

of Wronskians vanishes: 







cos 




el = 


2^-"e 


. sin 


fe = 


4>ee 


sin - 




n = 


277-"« 


. cos 



= W^e(e°,et(0)) + iye(/°,/t(0)) 



= e 



(e\ - cos-i/'eC^) - eg . (e^ - cosV'es") 



sinf/'e sinV', 



+ (/.^ - cos - /."-^ (/^ - COS 

^« W /.O „l^ , V'/ rr, / .0 .In 



sin^e sin^/ 
= -^(0.. f-^ cos COS + ^ sin is!^ sin ^\ 

_ ^,^1 f ^ cos i^^:^ sin - ^ sin i^^:^ cos ^\ 

+ ^2(^0 „ ^) eos - cos (2!£^ sin ^) | 

This is an ordinary differential equation for (p which is coupled to the (integral) 
equations for calculating the /3's as functions of (p. If it is non-singular (i.e., the 
coefficient function of (pge does not vanish for any 6*) then there is a solution 
at least locally. But the non-existence of the inf described for open curves above 
will also affect closed curves and global solutions may actually not exist. However, 
for closed curves that do not double back on themselves too much, as we will see, 
geodesies do seem to usually exist. 



4.5. An Example. Geodesies in the sphere are great circles, which go all the 
way around the sphere and are always closed geodesies. In the case of the Grass- 
mannian, using the Jordan angle basis, the geodesic can be continued indefinitely 
using formula 5 above. In fact, it will be a closed geodesic if the Jordan angles 
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ipe^ipf are commensurable. It is interesting to look at an example to see what 
sort of immersed curves arise, for example, at the antipodes to the point repre- 
senting the unit circle. To do this, we take (P{0) = e'^/27r to be the circle of unit 
length, giving the orthonormal basis e° = cos{9/2)/^, p = sin(^/2)/y^. We 
want e^, to lie in a direction horizontal with respect to these and the simplest 
choice satisfying the Wronskian condition is: 

/cos(20) 



The result is shown in figure 3. 



i sm{26) 







Figure 3. A great circle geodesic on Bod- The geodesic begins at the 
circle at the top left, runs from left to right, then to the second row and 
finally the third. It leaves Bod twice: at the top right and bottom left, in 
both of which the singularity of figure 1 occurs in 2 places. The index of the 
curve changes from +1 to —3 in the middle row. See text. 



5. Sectional curvature 

We compute, in this section, the sectional curvature of i3i^od/(sim) (i.e., trans- 
lations, rotations, scaling). Wc first compute the sectional curvature on the Grass- 
mannian which is non-negative (but vanishes on many planes) and conclude from 
O'Neill's formula [14] that the sectional curvature on Bj/(sim) is non-negative But 
since the O'Neill correction term is difficult to compute in this setting wc also do 
it in a more explicit way, computing first the curvature on the Stiefel manifold 
by Gauss' equation, then carrying it over to Imm /(transl). Since this is an open 
subset in a Frcchct space, the O'Neill correction term can be computed more easily 
on Imm / (transl) and so we finally get a more explicit formula for the sectional 
curvature on Bj/ (transl). 

5.1. Sectional curvature on Gr(2, F). Let W e Gr(2,y) be a fixed 2-plane 
which we identify again with M^. Let r} : V ^ V he the isomorphism which 
equals —lonW and 1 on W-^ satisfying r} = ri~^. Then the Grassmanian is 
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the symmetric space 0{V)/{0{W) x 0{W^)) with the involutive automorphism 
a : 0{V) 0{V) given by a{U) = rj.Ur]. For the Lie algebra in the V = W®W-^- 
decomposition we have 



-1 
1 



-y 
u 



-1 
1 



y u 



Here x G L{W,W),y G L{W,W^). The fixed point group is 0{VY 
0{W^). The reductive decomposition g = t + p is given by 



0{W) X 



X -y^ 

y U 



X 

u 



, X e so(2) 



,yeLiW,W^) 



Let TT : 0{V) ^ 0{V)/{0{W) x 0{W-^)) = Gr(2, V) be the quotient projection. 
Then T^tt : p To Gr is an isomorphism, and the 0(y)-invariant Riemannian 
metric on Gr(2, V) is given by 



■tr 





yi 



T 

-yi 





2/2 



T 

-2/2 





-2/12/J 



i trw'lyf y2) + h trw± {yivl) 



= ^^w{yly2) = {yi{ei),y2{ei))w^ + (2/1(62), 2/2 (e2))H'-L 

for Yi, Y2 e p, where ei, 62 is an orthonormal base of W. By the general theory of 
symmetric spaces [6] , the curvature is given by 

Rf'^inn.Y,, T,n.Y2)Ten.Y, = nn.[[Yr,Y2],Y,] 





2/1 

2/1 



-yl 






2/2 



-yl 




-2/fy2 + 2/^ 2/1 




-2/l2/^ + 2/22/f 



T 
-2/1 









"0"')] 


V2/2 









-2/f2/2 + yl 2/1 









-2/l2/^ + 2/22/? 







-2/f 
2/1 

22/^2/22/? - 2/^2/12/? - 2/f 2/i2/J 



v,-22/i2/?2/i + 2/22/^2/1 + 2/i2/f2/2 
For the sectional curvature we have (where we assume that Yi, Y2 is orthonormal): 

-B{Y2, [[Yi,Y2],Yi]) = trw{yly2ylyi + 2/ J 2/1 2/^2/2 - '^y'^yiylyi) 



, Gr(2,y) 
%pan(Yi,r2) 



= ^ trvK ((2/^2/1 - 2/^2/2)^(2/^2/1 - 2/^2/2)) 



2^2 

+ 5 trjyi ((^22/? - yiyl f{y2yl - yiyl)) 
= hhlyi - yiy2\\L^(w^w) + \\\y2yl -yivlfL^^w^^w^) ^ O- 

where L'^ stands for the space of Hilbert-Schmidt operators. Note that there are 
many orthonormal pairs Yi , Y2 on which sectional curvature vanishes and that its 
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maximum value 2 is attained when yi are isometries and 2/2 = Jyi where J is 
rotation through angle 7r/2 in the image plane of t/i. 



5.2. Sectional curvature on Imm/(sim). The curvature formula can be rewrit- 
ten by 'lowering the indices' which will make it much easier to express it in terms 

of the immersion c. Fix an ortho normal basis e, f oiW and let Sck — yfe(e), Sfk 
i. J . For x,y £ W-^, we use the notation xt\y = x®y — y®x & W-^ ® W-^. 
Then: 

CS^Iy.) = (('^'^1' '5/2) - (<5e2, 5h)f + ijlJei A Se^ + Sf, A Shf. 

To check this, note that yjyi — y'[y2 is given by a skew-symmetric 2x2 matrix 
whose off diagonal entry is just {Sei,Sf2) — {Se2,Sfi) and this identifies the first 
terms in the two formulas for k. On the other hand, y2yf is given by a matrix of 
rank 2 on the infinite dimensional space W-^. In view of C L{W-^, W-^) 

it is the 2-tensor 6ei ® Se2 + Sfi (g) 6/2- Skew-symmetrizing, we identify the second 
terms in the two expressions for k. 

Going over to the immersion c, the tangent vector Sck + iSfk to Gr becomes 
the tangent vector hk = 5c = J {5ek + iSfk){e + if)d9 to Imm/(sim). To express 
the first term in the curvature, we have: 

Proposition. {5e\,5f2) — {5e2,5f\) = f^det{Dshi,Dsh2)ds. 



Proof. 

n ih ^ (e + i/).(fcfc + i(5/fc) 

^^y^-k) = g2 + y2 

hence 

det{Dshi,Dsh2) = lm{Dshi, Dsh2) = ^^^^ 2 , f2 

hence 

/ det{Dshi,Dsh2)ds= [ {SeiSf2 - 5e25fi)d9. 
Jc Js^ 



□ 



The second term is not quite so compact: because it is a norm on W-^ (S> W-^, 
it requires double integrals over C x C, not just a simple integral over C. We use 
the notation as above c{6) = r(^)e*"(^'. Then we have: 

Proposition. 

\\Sei A Se2 + 5 fx A 6/2]^ = terml -|- term2 

terml -yy ^ \ -{D,h2{x),DM{y)) ) '^'^'^''^'^y' 

CxC 
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ff l-cos{a{x)-a{y)) f det(I);Z;i(-)' f ^Z^f )) ^ yds{x)dsiy) 



CxC 



Proof. Using r and a, we have: ^/re ^"^^Dghk = Sck + iSfk, hence: 

Skew-symmetrizing in the 2 vectors hi, /i2, we get: 

^r{x)r{y)Re {e""''^^"""'" (^D,hi{x)D,h2{y) - DJ^D,hi{y)j } = 

5ei{x)de2{y) - Se2{x)dei{y) + 6fi{x)5f2{y) - 5f2{x)6fi{y). 

Squaring and integrating over 5*^ x 5^, the riglit hand becomes \\5ei A 5e2 + Sfi A 
(5/2|p. On the left, first write Re(e*("(^)-"(''»/2(-. . . )) as the sum of cos((a(a;) - 
a(y))/2)Re(- • • ) and — sin((Q:(x) — Q!(y))/2)Im(- • • ). Then when we square and 
integrate, the cross term drops out because it is odd when x, y are reversed. □ 

We therefore obtain the expression of the curvature in Imm/(sim): 

CxC 

rr i-cosiaix)-aiy)) f detp^^,(x), I. M.)) ^ Y dsix)dsiy) 
77 2 V -det{Dsh2{x),Dshi{y)) J ^ y 

CxC 

A major consequence of the calculation for the curvature on the Grassmannian 

is: 

5.3. Theorem. The sectional curvature on Bi/ (sim) is > 0. 

Proof. We apply O'Neill's formula [14] to the Riemannian submersion 
TT : Gr° ^ Gr° /U{V) ^ Bi/ Diff+(6'^) 

where X^°'' is a horizontal vector field projecting to a vector field X at Tr{W); 
similarly for Y. The horizontal and vertical projections exist and are pseudo 
differential operators, see 5.6. □ 
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5.4. Sectional curvature on St{2,V). The Stiefel manifold is not a symmetric 
space (as the Grassmannian); it is a homogeneous Riemannian manifold. This can 
be used to compute its sectional curvature. But the following procedure is simpler: 

For (e, /) e we consider the functions 

QiieJ) = l\\ef, Q2{e,f) = \\\ff and ^3(6,/) = ^(e, /)• 

Then St (2, V) is the codimension 3 submanifold of V"^ defined by the equations 
Qi=Qi = 1/2, Qz = 0. 

The metric on St{2,V) is induced by the metric on V^. If ^1 = (5ei,(5/i) 
and ^2 = (<5e2,<5/2) are tangent vectors at a point in V"^, we have (^1, ^2) = 
((5ei , de^) + {5f\ , 6/2)- For a function on its gradient gradt^ of (if it 
exists) is given by (grad<p(u), ^) = d(fi{v){^) = Dy^^(p. The following are the 
gradients of Qf. 

gradQi = (e,0), gradQ2 = (0, /), gradQs = -^(/, e), 

and these form an orthonormal basis of the normal bundle Nor(St) of St{2,V). 

Let ^1,^2 be two normal unit vectors tangent to St{2,V) at a point (e,/). Since 
is flat the sectional curvature of St (2, V) is given by the Gauss formula [4]: 

Cn(L.) = (^(^i'^i)'^(?2,e2)> - (^(a,6),^(ei,6)) 

where S denotes the second fundamental form of St(2, V) in V"^. Moreover, when 
a manifold is given as the zeros of functions Fk in a flat ambient space whose 
gradients are orthonormal, the second fundamental form is given by: 

S{X, Y) = J2 Hf, {X, Y) ■ gvBAFk 

k 

where H is the Hessian of second derivatives. Given ^1,^2 G T'(e,/) St with = 
{Sei,6fi), we have: 

HqA^u^2) = {Se^,5e2), ifQ,(a,6) = ('^/i,^/2), 
HqA^u^2) = ^{{Seu6f2) + {662, Sh)) 

so that 

'5(6,6) = -(^ei,^e2)grad(5i - (^/i, (5/2) gradQ2 

- 7l(('5/i,^e2) + {5euSf2)) gradQa- 

Finally, the sectional curvature o/St(2, y) for a normal pair of unit vectors ^,r] 
in Tf St(2, V) is given by: 

CSL.) = Il^ei|ni^e2||^ + \\5f4'\\Sf2f + 2{6eu5f{){5e2,5f2) 

- {6ei,6e2f - {6fi,5f2f - i((5ei,5/2) + {Sfi,6e2)f 
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(10) = MSei ® 6e2 - 5e2 O 6ei + Sfi ® 6/2 - 5f2 ® 5fi 



|2 



U{6e,,6f2)-{Se2,6h)) 



2 



2 

Comparing this with the curvature for the Grasmannian, we see that the O'Neill 
factor in this case is ^{{661,6/2) + (<5/i, <5e2))^. Moreover, we can write for the 
curvature of the isometric Imm/ (transl, seal) 

, Imm/(transl.scal) if / j./t-it i \ i 

CxC 

+ 1 // I - - ^(^)) . ( det(D I^.;^2y )- yrf,(^)rf,(^) 

2 V -det(L',/i2(a;),Li,/ii(y)) ^ v ; vy; 



CxC 



5.5. Sectional curvature on the unsealed Stiefel manifold. Using the basic 
mapping the manifold Imm /(transl) can be identified with the unsealed Stiefel 
manifold which we view as the following submanifold of V'^. We do not introduce 

a systematic notation for it. 

(12) M = {{eJ)&V^\ {(0, 0)}, ||ef = and (e , /) = 0} 

equipped with the metric 

IK''e,<);)||(^j) - I ||g||2^ ||y||2 • 

Consider the diffeomorphism * : IR+ x St(2, V) ^ M defined by 

^{^,{e,f)) = {^rLe,^^Lf)=■. (e,/). 
For ^(5e, 5f) e T^^j) St we have 

%e,/)*.(A, = (^e + Vl.5e, ^/ + VL5f) =: {6e, 6f) 

Thus, ^ is an isometry if M+ x St(2, V) is equipped with the metric 

ll(^'Olll(e,/) = ^ + ll^e||2 + ||5/||2 

so that M is isometric to the Riemannian product of M+ and St(2,y), taking 

||A||^ — \\\/{^/2.£) for the metric on R+. This implies that the curvature tensor 
on M is the sum of the tensors on R+ (which vanishes) and St(2, y). Thus, if 
= T(£j-)7r.(A, ^j), i = 1,2 with (Ci,Ci) orthonormal, 

, M ^ -(i?M(a,6)6,6) ^ -(ii;St(6,6)6,6) 

span(a,6) ||^-i||2||^-2||2 _ (^-^,^-,)2 ||^-,||2||^-2||2 _ ^^^^^^^2 
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span({i,{2) 



Note that we have the relations: 

(13) 5e = -^e + yfl.5e 

(14) 



5.6. O'Neill's formula. For Riemannian submersions, O'Neill formula [14] states 

that the sectional curvature, in the plane generated by two horizontal vectors, is 
given by the curvature computed on the space "above" plus a positive correction 
term given by (3/4) times the squared norm of the vertical projection of the Lie 

bracket of any horizontal extensions of the two vectors. We now proceed to the 
computation of this correction for the submersion from Imm /(sim) to _Bi/(sim). 

Because of the simplicity of local charts there, it will be easier to start from 
Imm /(transl). Let c G Imm with J^j cds = 0. We first compute the vertical 

projection of a vector h G Tc Imm /(transl) for the submersion Imm /(transl) — > 
i?j/(sim). Vectors in the vertical space at c take the form 

h = bv + iac + (3c, 

each generator corresponding (in this order) to the action of diffeomorphisms, 
rotation and scaling {b is a function and a,/3 arc constants). Denoting hJ the 
vertical projection of and using the fact that Gc{h,h) = Gc{hJ,h) for any 
vertical h, we easily obtain the fact that 

h'' = bv + iac + f3c, 

with 

L^b + an = v ■ Lh, 
{bn) + a = {Dgh ■ n), 
13 = {Dsh-v), 

where we have used the following notation: Lh = —D^h, L^b = —D^b + K^b and, 
as before 

jF) = j J Fds. 

From this, we obtain the fact that b must satisfy 
(15) L^b - Jbi^jK ^vLh- {Dsh ■ n)K. 

The operator is of order two, unbounded, selfadjoint, and positive on {/ e 
L'^{S^,ds) : Jfds = 0} thus it is invertible on {f G C°° {S^ ,R) : J f ds = 0} by 
an index argument as given in lemma [11, 4.5]. The operator in the left-hand 
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side of (15) is also invertible under the condition that c is not a circle, with an 
inverse given by 



(16) (i^)- V = (i^)- V + -^^^^=^-^ii7r'n. 

1 — {k[L ' )^k) 

This is well defined unless k = constant. Indeed, letting / = {L^)~^k, we have 
—fD^f + K^f'^ — kJ which implies {nf) > (k;^/^). By Schwartz inequality we 
have < iWP))^^'^ which ensures (nf) < 1. Equality requires {fD^f) = or 
/ — constant, which in turn implies that k = constant and that c is a circle. We 
note for future use that (Z^)^^k — {L^)^^k. 

We hereafter assume that c has length 1, is parametrized with its arc-length 
divided by 2n, and that it is diflferent from the unit circle (which is a singular 
point in Bi/(sim)). We can therefore write 

(17) = (^(Z^)-V(/i)w-»(«(^'^)"V(/i))c^ +i{Dsh-v)c+ {Dsh-n)c 

with ip{h) = V ■ Lh— {Dgh ■ n)K. 

The right-hand term in (17) is the sum of three orthogonal terms, the last two 
forming the vertical projection for the submersion Imm /(transl) Imm /(sim). 
Applying O'Neill's formula two times to this submersion and to Imm/(sim) — » 
Bi/(sim), we see that the correcting term for the sectional curvature on Bi/(sim) 
relative to the curvature on Imm /(sim), in the direction of the horizontal vectors 
hi and /i2 is 



p{hi,h2)c = ^ 



(L^)-V([/ii, h2]c)v - linii-^y^^iihiMc 

h\,h2 being horizontal extensions of h\ and /i2- From the identity 

= j {bLb + K%'^n - {Kb)Kb)ds 
= j b{L^b)ds 

we can write 

p{hi,h2)c = \ j i!{[huh2]c){L')-^i^{[hiMc)ds. 
We now proceed to the computation of the Lie bracket: 



5.7. Proposition, ■^([/ij^, /i^]c) = Ws{Dsh\ ■ n,Dsh2 ■ n) — {det{Dshi, Dsh2))K 
where Ws{h,k) = hD^k — kDgh is the Wronskian with respect to the arc-length 
parameter. 
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Proof. We take ft.1,/12 e {/ G C°°(5\M2) : J f ds ^ 0} which are horizontal at 
c, consider them as constant vector fields on Imm /(transl) and take, as horizontal 
extensions, their horizontal projections 7 h^il), ^2 (7)- Then we compute the 
Lie bracket evaluated at 7: 

[hi, hi] = D,^i,Xh)-DcMhih) 

= -D^.h.hJi-fj+D^^h.hiij) 
since hj + = hi is constant, for i = 1, 2. We have 

+ i{D,_^hi ■ n^) j + (Ds^hi ■ v^) "f 



with xl)^(h\) = v-y ■ Ljhi — (^Ds^h\ ■Vj)^Kj. We have added subscripts 7 to 
quantities that depend on the curve, with Dg^ holding for the derivative with 
respect to the 7 arc- length (we still use no subscript for 7 = 0). Note that 
(^D s^hi ■ n^') ^ = £^{Dshi-n^/) and (^Dg^.hi ■ v^'^^ = i^{Dshi ■ v^) which is a 
first simplification. Also, since we assume that hi is horizontal at c, we have 
{Dghi ■ n) = {Dghi ■v)=0 and v ■ Lh\ = 0, which imply V'(/ii) = 0. 

We therefore have (to simplify, wc temporarily use the notation /' = Dgf) 
(18) D,,hX = (^{L'^)-^DeM^^{hi)v-i{K{L^)-'D,,hMhi))c^ 
+ i{h[ ■ Dc,h2'nj)c+ {h[ ■ Dc,h2V'y)c. 

Since we have Dch^v-y = {h'2 ■ n)n and Dch^n-^ — —{h'2 ■ 'n)v we immediately 
obtain the expression of the last two terms in (18), which are 



(19) -i{{h'x ■ v){h'2 ■ n))c+{{h[ ■ n){h'2 ■ n))c. 

We now focus on the variation of (f)^. Wc need to compute 

Dc,h2'^J-f{hi) = DcMiv-f ■ L^hi) - {{h[ ■ v){h2- n))K. 
If /i is a constant vector field, we have 

Ds,h = h'\W\\-' 

and 

This implies 



L,h = -ih'\h'\\-'ytf'\\-'. 



D^h^L-yhi = -/i'/£'c,fcj|7'll"' - {h[Dc,h2\\D,-f\\-^y 
= 2h'({h'2 ■ v) + h[{h'2 ■ v)' . 
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Therefore 

De,h,{L^h ■ v^) = (/I'l • v){h'2 ■ v)' - {h'l ■ n){h'^ ■ n). 

Using 

h'l = {{h\ ■ v)v + {h[ ■ n)n)' 

= ((/I'l • v)' - K{h[ ■ n))v + {{h[ ■ n)' + K{h[ ■ v))n 
and the fact that h'l ■ v = h'2 ■ v = Q, we can write 

Dc,h2{Ljhi ■ Vj) = -(/12 • n)(/ii • n)' 

which yields 



D^,h,Mf^i) = -{h'2 ■ n){h'^ ■ n)' - ((/I'l • v){h'2 ■ n))K. 

By symmetry 



DcM'^ii^i) - Dc,hii'i{h2) = VFs(/ii • n, h'2 ■ n) - (det(/ii, /iy)K;, 
where Ws{ipi,ip2) = ^pw'2 - 

Combining this with (19), we get 



\ht,h^\c = {L^y^{Ws(h:^-n,h!2-n)- (Aci{h!^,h!2))K)v 



K{p)-^{Ws{h'^ ■ n,h'2 ■ n) - (det(/i'i, /i'2))k) )c - i(det(/i;, ft'2))c 

so that 



^jilh^h:^],) = •n,fc^-n)-(det(/i'i,/iy)K. 
(We have used the fact that ip{bv + iac) = L^b.) □ 

We therefore obtain the formula 

(20) p{hx ,h2)c = ^J(Ws{D,hi-n,D,h2-n)- (det(£>,/ii, £>,/i2))k) • 

• iL^y^(Ws{Dshi-n,Dsh2-n) - {det{Dshi, Dsh2))Kys 
with (I/^)"^ given by (16). Finally, assuming that hi and /12 are orthogonal, 

,Bi/(sim) _ , Imm/(sim) , u \ 

wliere fc^^f^™^ gi^^n in (9). 

A similar (and simpler) computation provides the correcting term for the space 
Bi/(transl, scale). In this case, the rotation part of the vertical space disappears, 
and the two remaining components (parametrization and scale) are orthogonal. 

The result is 

(21) p(/ii,/i2)c = 7 / WsiD,hi-n,D,h2-n){L'^r^WsiD,hi-n,D,h2-n)ds. 
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5.8. An upper bound for fc^'al(™/i2)' Here we derive an explicit upper bound 
for k^'^^T^\ N at a fixed curve c e Bj/fsim) and a fixed tangent vector /lo. This 
will show that geodesies (such as the one in the hi direction) have at least a 
small interval before they meet another geodesic. The size of this interval can be 
controlled, as we will see, by an upper bound that involves the supremum norm 
of the first two derivatives of hi. 

We assume that c has length 27r. Since Imm/(sim) is isometric to Gr{2,V), 
its sectional curvature is not larger than 2 as already remarked. We estimate the 

terms in p(ft.i,/i2)c = /i2)(£^)"V(^i, ^2) > where ^'(^1,^2) = WsiDghi ■ 



n,Dgh2 ■ n) — {dct{Dshi, Dsh2))K. For a fixed ft,2, '>P{hi,h2) is function of hi be- 
longing to H~^{c). We estimate /i2)||c,-i and then /i2)c by estimating 
the norm of the operator which maps i?~^(c)to H^{c). 

Iff€ H\c), ||/||e,-i < ll/lico and ||/'||e,-i < ||/||c,o. Therefore, 
\\Ws{D,hi-n,D,h2-n)\\_i = \\{D,hi-n)D,{D,h2-n)-D,{D,hi-n){DM-n)\\-i 
< {\\Dsh2 ■ n||c,cx> + \\Ds{Dsh2 ■ n)\\c,ao) \\Dshi ■ n||c,o- 
Since hi has norm 1, Hfs/ii • n||c,o and ||Z)s/ii • wUcO are < \/27r. 



{det{DM,D,h2)) < {\D,hi ■ n\ ■ \D,h2 ■ v\) + {\DM ■ v\ ■ \Dsh2 ■ n\) 

< {\\D,hi ■ n||e,o • \\Dsh2 ■ v\\e,o + \\D,hi ■ v\\e,o ■ \\D,h2 ■ n\l,o) < 2. 
This results in 

h2)\\c,-i < V2^ (^\\Dsh2 ■ n||e,oo + \\D,{D,h2 ■ n)||e,oo + 2y^^ . 

Now 



1 - - 1 _ (k(lT)-i^^ 



since (Vj(LT)-i«;) < . {^{L^)-^^). 

5.9. Proposition. If tpeH~^{c) then 



27r 



-1- 



Proof. Let = + 1. Let L'^ f = LJo = ^- Then, fJoeH\c) and 
||/o||c,i = IIV'llc-i- Let 5 = / -/o so that L'^g = (1-k^)/o. The eigenvalues of L"^ 
are positive and bounded from below by 1/2, see [3]. Therefore, j|g||^ g < 2{g, L^g) 
where {g,L^g) = J gL^gds. We also have Hg'ljco ^ {9i^^ 9) Hence, 



< ^9,L'9) < 3/(1 - i^'')9fods < 3||1 - «2||c,oo||5l|c,i||/o| 



c,l- 
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Therefore, \\g\\e,i < 3|| 1 - and ||/||c,i < (l + 3|| 1 - ||^_^) 

Finally, 

MI^PW < ^ll^llc-i • ll(i^)- Vllci < ^ (1 + 3||1 - kIcoo) □ 



Putting all the estimates together we get, for orthonormal hi, /i.2 as always, 

; 
(22) 



" - ''^spaii(/ii,ft2) - 



< 2 



3 (1 + 3||1 - k'^Wcoo) (^\\Dsh2 ■ n\\c,oo + \\{Dsh2 
4(l-(/.(iT)-i«)) 



6. Numerical procedure and experiments 

The distance -Dop.dif given in 4.1 can be computed in a very short time by 
dynamic programming, using a slightly modified procedure from the one described 
in [16]. Here is a sketch of how it works. 

Let i^(a",a^) = max(0, cos((q!° — a^)/2)), and assume that the curves are 
discretized over intervals [6^ {k) , 0^ {k + 1)) , fc = 0, — 1, i = 0, 1, so that the angles 
have constant values, a* (A;) on these intervals. The problem is then equivalent to 
maximizing 

,-min{e°{k+l),<p-'^{e^{l+l))) 
" ./max(eO(fc), 0-1(^1(0)) 



1.1 Jm; 



with Fki = -F(a°(fc), a^(l)). Because the integral of is maximal for linear </>, 
we must in fact maximize 

Y^Pklsj (ma^ieo {k),e^m - min(^0(fc + l),^i(/ + 1))))^ 



(max (9°{k), 9^1))) - mm{e%k + 1), 6^1 + 1)))) ^ 

with the notation 6i°(fc) = (/)(6'°(A;)) and 9^(1) = (f)-^{e\l)). The method now 
essentially implements a coupled linear programming procedure over the values of 
^0 and Oi. See [18, 16] for more details. This procedure is very fast, and one still 
obtains an efficient procedure by combining it with an exhaustive search for an 
optimal rotation. 
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For closed curves, we can furthermore optimize the result with respect to the 
offset 0(0) e , for the diffeomorphism. Doing so provides the value of 

£';p,diff,rot(c°,ci) = infarccosV(C-(.^))2 + (S_(.^))2. 

where the notation D' is here to remember that the minimization is over </> e 
C°°'+{S^) and not C°°'+([0, 27r]). 

This combination of the almost instantaneous dynamic programming method 
and of an exhaustive search over two parameters provides a feasible elastic match- 
ing method for closed curves. But this does not provide the geodesic distance 
over Bi/ (sim), since wc worked with great circles instead of the Neretin geodesies. 
There are two consequences for this: first, the obtained distance is only a lower 
bound of the distance on B^/ (sim), and second, since the closedness constraint is 
not incliidc;(l, the curves generally become open during the evolution (as shown in 
the experiments). 

However, the optimal diffeomorphism which has been obtained by this approach 
can be used to reparametrize the curve c^, and we can compute the geodesic 
between c° o 0* and in Imm/(sim) using Neretin geodesies, which forms, this 
time, an evolution of closed curves. Its geodesic length now obviously provides an 
upper-bound for the geodesic distance on Bi / (sim). The numerical results that 
are presented in figures 4 to 8 compare the great circles and Neretin geodesies 
obtained using this method. Quite surprisingly, the differences between the lower 
and upper bounds in these examples are quite small. 




Figure 4. Curve evolution with and without the closedness 
constraint. Lower and upper bounds for the geodesic distance: 
0.443 and 0.444 
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Figure 5. Curve evolution with and without the closedness 
constraint. Lower and upper bounds for the geodesic distance: 
0.462 and 0.464 




Figure 6. Curve evolution with and without the closedness 
constraint. Lower and upper bounds for the geodesic distance: 
0.433 and 0.439 




Figure 7. Curve evolution with and without the closedness 
constraint. Lower and upper bounds for the geodesic distance: 
0.498 and 0.532 



7. Appendix: The geodesic equation on G!™m,scai,i, 



7.1. The geodesic equation. We use the method of [11] for the space {c G 
Imnic : c(l) = 0} which is an open subset in a Frechet space, with tangent space 
{h e C°°(5^,C) : /(I) = 0}. We shall use the following conventions and results 
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Figure 8. Curve evolution with and without the closedness 
constraint. Lower and upper bounds for the geodesic distance: 
0.513 and 0.528 



from [11]: 

= I?s,c = j^, ds=\ce\d9, Dc,rJ-{c) = J {Dsm,Vc)ds = - J Kc{m,nc) ds, 

Dc,m{Ds) = -{Ds'm,Vc)Ds, Dc,mids) = {Dsm,Vc)ds. 
Then the derivative of the metric at c in direction m is: 

Dc,mGc{h, k) = ^ j K,c{m, ric) ds ■ J {-D'^h, k) ds 

+ j J{D,m,Vc){D% k)ds+j j (D,{^{Dsm,Vc)Dsh^,k'^ ds 
-jj {Dlh,k){D,m,Vc)ds 
= J2 J '^c{m, He) ds ■ J {-D^h, k) ds 
-jj (^-Dlm,D-^[{D,h,D,k)v^)ds 

According to [11, 2.1] we should rewrite this as 

fc) = Ge(Jr?(m, /i), fc) = Ge(m, fc)) , 
and thus we find the two versions K and H of the G-gradient of c i— > Gc{h, k): 

Kc{m, ^ Y J '^"^^^ '"'"^ ds ■ h - D^^ {{Dsm, Vc)Dsh) 

Hc{h,k) = -^D^^K^Uc) ■ J {-Dlh,k)ds-D-^({D,h,D,k)vo) 

which gives us the geodesic equation by [11, 2.4]: 
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(23) ~ T J "'"'^ ds-ct- D^^ {{DsCt, Vc)DsCt) 

7.2. Theorem. For each k > 3/2 the geodesic equation derived in 7.1 has unique 
local solutions in the Sobolev space of H'' -immersions. The solutions depend C°° 
on t and on the initial conditions c(0, . ) and ct(0, . ). The domain of existence 
(in t) is uniform in k and thus this also holds in Imm, := {c e Imm(S'^,IR^) : 
c(l)=0}. 



Proof. The proof is very similar to the one of [11, 4.3]. We denote by * any space of 
based loops (c(l) = 0). We consider the geodesic equation as the flow equation of a 
smooth (C^) vector field on the H^-opcn set U'^ x M^) in the Sobolev space 

H^{S^,R^)xH^{S\R^) where C/*^ = {c e : \ce\ > 0} C H'' is F^.open. To see 
that this works we will use the following facts: By the Sobolev inequality we have 
a bounded linear embedding H^{S^,M.^) C Cr(S\ M^) if fc > to+ i. The Sobolev 
space H^{S^,M.) is a Banach algebra under pointwise multiplication if A: > ^. For 
any fixed smooth mapping / the mapping u i— > / o u is smooth H'^ if fc > 0. 

We write Ds^c '■= Dg just for the remainder of this proof to stress the dependence on 
c. The mapping (c, u) ^ —D^ c'^ is smooth U x ^k-2n jg ^ bibounded 

linear isomorphism — > iJ^-^n j^j. f^xed c. This can be seen as follows (compare 
with [11, 4.5]): It is true if c is parametrized by arclcngth (look at it in the space 
of Fourier coefficients). The index is invariant under continuous deformations of 
elliptic operators of fixed degree, so the index of —D^ is zero in general. But —D^ 
is self-adjoint positive, so it is injcctive with vanishing index, thus surjective. By 
the open mapping theorem it is then bibounded. Moreover {c,w) ^ {—D^)~^{w) 
is smooth U'^ x iJ^-2n _^ jjk ^j^g inverse function theorem on Banach spaces). 
The mapping (c, /) ^ DJ = j^^e/ is smooth x iJf 3 /7 x iJ™ ^ F™"! 
for k > m, and is linear in /. We have v = Dg.cC and n = iDg cC The mapping 
c I— > k(c) is smooth on the ff^-open set {c : \ce\ > 0} C iJ^ into H^~^ . Keeping 
all this in mind we now write the geodesic equation (23) as follows: 



Ct = u=: Xi{c,u) 

Ut = -D~'^(^^\\u\t=o\\G-l^c.nc + ^Ds^ci\Ds,cuf.Vc) + Ds,c{.{Ds,cU,Vc).Ds,cU) 

-{Yj{^,Dl,c)ds) 



t=0 



=: X2(c, u) 



Here we used that along any geodesic the norm jlctHc and the scaling momentum 
— J- J{ct,D^^^c) ds = dt\og£{c) are both constant in t. Now a term by term 
investigation shows that the expression in the brackets is smooth U'^ x H'' — > H''~'^ 
since k — 2 > |. The operator —D~1 then takes it smoothly back to H'^. So the 
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vector field X = {Xi,X2) is smooth on U'' x H''. Tlius tiie flow Fl'' exists on H'' 
and is smooth in t and the initial conditions for fixed k. 

Now we consider smooth initial conditions cq = c(0, . ) and uq = Ct{0, . ) = 

7/(0, . ) in C°°(S'\K2). Suppose the trajectory Flf(co,uo) of X through these 
intial conditions in H'' maximally exists for t € {—ctk^bk), and the trajectory 
Flj"'"^(co, Wo) in -ff*^"*"^ maximally exists for t G {—ak+i,bk+i) with bk+i < bk- By 
uniqueness wc have Fl'^'^^ (cq , uq) = Flf(co,uo) for t £ {—ak+i.bk+i)- We now 
apply de to the equation ut = X2{c, u) = —Dj1{ . . . ), note that the commutator 
[de, —Dj^] is a pseudo differential operator of order —2 again, and write w = deu. 
We obtain wt = dgUt = —D~'^dg{ . . . ) + [dg, —D~'^]{ . . . ) + const. w. In the term 
dg{. . .) wc consider now only the terms dgU and rename them dgW. Then we get 
an equation wt{t,0) = X2{t,w{t,0)) which is inhomogeneous bounded linear in 
w e H'^ with coefficients bounded linear operators on which are C°° functions 
of c, u € H'^. These we already know on the intervall {—ak,bk). This equation 
therefore has a solution w{t, . ) for all t for which the coefficients exists, thus for 
all t e (afc,6fc). The limit limtyb^_^^'w{t, . ) exists in H'^ and by continuity it 
equals dgu in H'^ at t = bk+i- Thus the iJ'^+^-flow was not maximal and can be 
continued. So {—ak+i,bk+i) = {ak,bk)- We can iterate this and conclude that the 
flow of X exists in n„>fe H"' = C°°. □ 
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